arXiv:1501.07384vl [physics.hist-ph] 29 Jan 2015 


The Cofounder of Quantum Field Theory: Pascual Jordan 


W. Dittrich 

Institut fiir Theoretische Physik 
Universitiit Tubingen 
Auf der Morgenstelle 14 
D-72076 Tiibingen 
Germany 

electronic address: qed.dittrich@uni-tuebingen.de 
Abstract 

A comparative study is undertaken that brings to light the two different methods of how to treat the 
many-body problem in quantum field theory. The two main researchers who published the first versions 
of how to quantize a many-body assembly were P. Jordan and P.A.M. Dirac. What they understood by 
the so-called “second quantization” will be the subject of the paper. We will argue that it is Jordan’s field 
operator approach that until now constitutes the basis of any work in quantum field theory. 
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INTRODUCTION 


We will begin with a detailed study of the continuous string. As a simple one-dimensinal model 
the string with fixed end points will enable us to introduce the quantization of a field as function 
of ordinary {x,t) space. Here, second quantization starts with an operator-field (f){x,t). Our goal 
is to provide a simplified model for computing the analogue of Einstein’s formula for the energy 
fluctuation of black body radiation. [l| This was the starting point of Jordan’s seminal contribution 
to quantum field theory (QFTh). His work was published in the last chapter of the famous “Three- 
Man Paper” (Drei-Manner-Arbeit) and for a long time criticized even by his coauthors. Jordan 
himself rated his calculation as “almost the most important contribution I ever made to quantum 
mechanics (Q.M.).” 


It is a pity that among the founders of Q.M. Pascual P. Jordan is the least known, although he 
contributed more than anybody else to the birth of QFTh. S.S. Schweber calls him “the unsung 
hero” among the creators of Q.M. ^ Just the opposite holds true for P.A.M. Dirac. Many experts 
in the field consider his paper of 1927 the beginning of quantum electrodynamics j^. While Jordan 
uses (like in classical field theories) the ordinary space-time coordinate {r^t) as arguments of the 
quantum field, Dirac starts with the Schrodinger function with multiparticle space coordinates as 
argument. The second quantization is then introduced after a series of canonical transformations 
instead of elevating the one-particle Schrodinger field to an operator field that has to satisfy certain 
equal-time communtation relations. However it was shown rather soon after Jordan’s and Dirac’s 
papers appeared that the two methods lead to identical results. Although the classical theories 
of particle and field picture are completely different as far as their mathematical and physical 
content is concerned, their quantized versions are totally identical. Q] This will be demonstrated 
in the final chapter in modern terms. It remains to say that Jordan’s original approach to second 
quantization, not Dirac’s, became the standard procedure among researchers and in textbooks to 
formulate QFTh., be it non-relativistic or relativistic, which is the accepted language of elementary 
particle physics. 


I. THE CLASSICAL VIBRATING STRING 


Let us consider a vibrating string with fixed endpoints at x = 0 and x = 1. Let its linear mass 
density be a and it tension be r, both of which are assumed to be constant. Then the equation of 
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motion for the vertical displacement u{x, t) at point x and time t is given by 


1 d‘^u{x,t) d'^u{x,t) 




dx"^ 


with ri(0, t) = 0 = u{l, t) 


( 1 . 1 ) 


c = denotes the velocity in the string. 

Clearly, equation (|l.ip is the wave equation in one dimension. (A derivation of equation (11.11) will 
be given later with the aid of the action principle.) 

To solve the wave equation let us try the Ansatz in terms of normal modes qn{t): 


u{x, t) = '^ Un{x)qn{t) = ^ Unix)cn cos(a;„t + (/>„ 
n=l n=l 

To satisfy the endpoint conditions u{0,t) = 0 = u{l,t) we find for Un{x) 

/ 2 \ mrx 

Wn(ic)=l^l sm—p, n = l,2,.... 

Un{x) is called the n-th wave function, which satisfies the equation 

d? 2\ / \ 

T-^ + aUj^jUn[x) =0, OJn = —C. 
The Un{x) form a complete set of orthonormal functions: 

OO 

Un{x)Un{x') = 5{x - x') 


n=l 


or 


/2 mrx j2 mrx' , 

Z^\It sin 7 sm = 6{x - x ), 


n=l 


i i \ i 


i 


L 

J Un{x)aUm{x)d: 


I 

2 f mrx rmrx 
0 

-a 

— / 2 — ^nm ■ 


dx 


Incidentally, the Green’s function of the string is given by 


G{x, x'\uj) = 


n=l 


with uJn = ^V 7 = or Vn = nfi- 


Un{x)Un{x') 

/ i2 , ,2 

UJr, — 


2 sin ^ sin ^ 


rrl < ^ 


crZ ^ 




G{x,x'\uj) given in (II.7p solves the differential equation 

d^ 


— ^ x'\uj) = 5{x — x') 


( 1 . 2 ) 


(1.3) 


(1.4) 


(1.5) 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 
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with G{Q,x'\u)) = O = G{l,x'\uj). 


The boundary conditions are obviously satisfied; acting with — on the Green’s func¬ 

tion solution (II.7p leads back to the completeness relation (I1.5I1 . 


Transforming back to space-time, 


G{x, x'] t — t') 



ei^it-nG{x,x'\uj) 


we find 


G{x, x'] t — t') = Un{x)— - — - '^-Un{x')Q{t — t') . 


(1.9) 


n=l 


The kinetic energy of the string is obtained from the kinetic energy for an element of the string, 
\{adx)u^, and integrating over the total length I. 




( 1 . 10 ) 


Using the Ansatz (11.21) we obtain with the aid of (11.61) 


T = 


2 X!/ QnQm I Un{x)aUm{x)d2 


n^m 

oo 


= 2^<in- 

71=1 

^ 2 

We also need the potential energy U = f ( ) dx or, after an integration by parts and 

fixed endpoints, u{0,t) = 0 = u{l,t) 


U = - 


t 

J U{x,t) 


d‘^u{x, t) 
dx"^ 


dx. 


( 1 . 11 ) 


Again using the Ansatz (II.2p and equation (11.41) . i.e., 


d 


-y^Un{x) = - 0J^Un{x) , 

dx^ T 


Y _^ 'Y f 

U = --T'^qn{t)qm{t){-i^n)- / Un{x)(JUm{x)dx 

n,m Q 

oo 


n=l 
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Therefore the Lagrangian is given by 


- OO 

L = T-U = = (1.12) 

n=l 

This result is quite remarkable. The normal coordinates reduce the Lagrangian to a discrete set of 
uncoupled harmonic oscillators. The equations of motion for the normal coordinates are 

gn + <^n9n = 0, U = 1, 2, . . . , OO . (1.13) 


The total energy is given hy E = T + U. If we substitute qn{t) = Cn cos {unt + then we obtain 

^ OO ^ OO 

r = - ^ qn{tf = (-Cn sin (uint + (j)n)f 

n=l n=l 

^ OO ^ OO 

^= 2^1 = 2 Z] • 

n=l n=l 

Hence the total energy E of the string is a constant of motion: 

^ OO 

E = T + U = -Y,u^yn. ( 1 . 14 ) 

n=l 

Again, it is quite remarkable that the total energy of the vibrating string is given by a sum of 
contributions from each of the normal modes. 

To formulate Hamilton’s principle for the string we need a Lagrangian, so that we can perform the 
variation: 

ITi2 = J dtL, 

t2 

SWu = G[ti] - G[t2], 


where 


L{t) = 


dxC I (p{x, t), 


d(p{x,t) dif{x,t) 


dt 


dx 


X, t 


(1.15) 


0 

(Notice that we replaced u{x,t) by ip{x,t).) 

The Lagrangian density C is itself a function of the “field” ip{x,t). The role of the generalized 
velocities is taken over by ^^^^.The dependence on is new; it represents an “interaction” 

between the infinitesimally separated elements of the system. Note that C has the dimension of 
energy per unit length. 

In general the generalized coordinate (p{x,t) is varied according to 


(p{x, t) ip{x, t) + d(fi(x, t) 
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with the spatial boundary condition 5ip{x = 0,t) = 0 = 5ip{x = Z,t). We will, however, also allow 
for the time variation at the endpoints: 


5ip{x,ti) ^ ^ 5ip{x,t2), for all x. 

Hamilton’s principle then takes the form 


where 


with 


*1 

6Wi2 ~ J 

*2 


i 

L{t) = J dxC 
0 


a /d(p\ T f dip 
~ 2 J ~ 2 

= T-V. 


, a,T = const. 


(1.16) 


Using 


with 


we obtain 


dC dC dC dC (dip\ dC (dip 

“ = « * + (w) + -pm' 


dt J 


’ 


d 


dC _ 

dip \ V 
dt 


d 


dx J 


dC _ _ 
dip \ ^ V 9a: 

dx 


After a few more steps we finally end up with the Lagrangian formulation: 


ti i 

dWi2 = f dt f dx 


t2 0 
- I O'— -TT - r 


a 


dip 

'm 


A 

dt 

d’^ip 92 ^X Q 

dt^ ^ dx^) ^ dt 


From here we can identify: (51Ui2 = 0) 
dip d'^ip 


^ 


dt"^ dx"^ 


= 0 


5t : 


A 

dt 


1 + -(^ 

2 \ dt 2 [dx 


J ^ [2\dtJ 2 \dxJ ] 


Z. (+ -(^ 

2 [dt 2 [dx 


equation of motion, i.e., one¬ 
dimensional string equation. 

energy conservation 


= 0 


(1.17) 


(1.18) 


6 



or 


|(T + V = £) = 0 . 


Surface term 


i 


G = j dx 
0 




with 


and 


T-r dip ^ 

n = cr—— , momentum density 

ot 


a {dip\^ T /dtp\^ 

+2 te) density 


L 

E = J dx£ . 

0 

The momentum density conjugate to ip{x,t) is defined as 

n(a:, t) = , ' (= cr— for the string) 

di^\ dt 


at J 

The Hamiltonian density is in general given by 


w|4,,n,gix,t)=n|^-£. 


For our one-dimensional string this yields with (11.161) and (I1.23P : 

dip dtp / <7 / dip\^ T ( dip\^\ a ( dip 


n = a 


dt dt 


2 \ dt / 2 \dx 

2 


2 \ dt 


+ 


dx J 




and therefore 


L 

H(t) = 


In the Hamiltonian formulation we obtain: 

1 i 


6 Wi 2 = J dt J dx — 7i5t] 


^ , du dn d dn 

+ -TiT- —+ 


dt dip dx 


.'dip dn\ , d^, 


(1.19) 


( 1 . 20 ) 


( 1 . 21 ) 


( 1 . 22 ) 


(1.23) 


(1.24) 


(1.25) 
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so that 5Wi2 = 0 implies 


Surface term: 


. _dn _ dn d &H 

^ ' dt dtp dx 


<5n 


5t 


dip 

'di 

dH 

~dt 


_ m 
~ M 

= 0 . 


G = 


L 

j dx {Jl{x^t)5p{x,t) —V.6t) , 
0 




(1.26) 


(1.27) 

(1.28) 


II. PASCUAL JORDAN’S FIELD QUANTIZATION OF THE STRING - THE ORIGIN 
OF QUANTUM FIELD THEORY 

After the introduction to classical Hamiltonian mechanics for fields, we can now begin to quantize 
our so far classical string field according to canonical quantization. For this reason Jordan replaced 
the classical fields p{x,t),Il{x,t) etc. by quantized operator fields which have to satisfy certain 
equal-time commutation relations. To derive these relations I will state without proof the following 
fundamental operator statement of the quantum action principle: Suppose we have an operator 
O = 0{p{t)) and we consider a variation corresponding to a change of the parametric variable t; 
then the infinitesimal change of the operator O is given by 

50 = -^[0,G]. (2.1) 

Let at a given time {5t = 0) the variation 6p{x,t) / 0; then the operator version of (11.271) is given 
by 

i 

G{t) = J dx'll{x',t)5p{x',t). 

0 

Now, using (j2.1[) with O = p{x,t), we have: 

5p{x,t) = 

8 


p{x,t), / dx'll{x',t)6p{x',t) 








or 


[(/^(x, t), n(a:', t)] =ih6{x — x'). 

Here is a list of Jordan’s equal-time commutation relations in the Heisenberg picture: 


( 2 . 2 ) 


Vp{x,t),ip{x' ,t) = fe,9i]=0, 

[n(x,t),n(x',t)] = 0, boPi] = o> 

[ip{x,t),Il{x' ,t)]= ih6{x - x'), [qi,pj] = ih6ij . (2.3) 


The discrete analogies are on the right. For the string field we found n(x,t) = 

In the following we want to employ our results obtained so far to compute the mean-square energy 


fluctuation in a narrow frequency interval + Aw) in a small segment of the string, i.e., in the 


region 0 < x < a, a <C b We follow mainly the paper by the author of ref. [^. Our former (p{x, t) is 


from now on considered a quantized field. Again we expand in slightly modified modern notation: 


oo . ^ y /2 / 2 

ip{x,t) = I —j g'„(t)ttn(x), where, as before , ttn(x) = ( — ) sinfe^x (2.4) 


n=l 


with 


qn{t) = , ujn = ckn = n . (2.5) 

At this point it would be useful to say a few words about the Hermitean particle operator N = a^a 
and the creation and annihilation operators and a of the excitations of the normal modes of 
the string. It will also help us to better understand Dirac’s definition of of bl and br in terms of 
the amplitude and phase (0^)- The great idea of Dirac was to make the Schrodinger wave 

function a multiparticle function of particle numbers N'^,'ip{N[,N 2 ,..., ■ ■ ■ ,t) and show how 

this function (not state) is changed when acted upon by the creation and destruction operator br 
and bl- Let us summarize some known facts from oscillator physics: 


So we obtain 


N = a^a = , [a, a'*^] = 1 

[a, iV] = aa^a — a^aa = (^aa^ — a = [a, a)]a = 


( 2 . 6 ) 


aN — Na = a , 
or aN{N + l)o , 
or a{N — 1) = Na . 


(2.7) 
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Likewise for the adjoint a^: 


a)N -No) = -a^ 
or No) = a'^(A^ + 1), 
or {N — l)a^ = a)N . 

Starting with aN = {N + l)a we obtain immediately aN'^ = {N + l)aA^ = {N + l)^a,... 
Hence, for every function that can be expanded in a power series, we have 


af{N) = f{N + l)a 


an example of which is 


where A is a parameter. The proof for this important result is simple: 


e^^ae = a + \[N, a] + a]] H- 

A \ -A 

— (1 — — - -{- —- — —- T * * * I n — c a ■ 
1! 2! 3! 


Therefore 


-\N -XN -A \ -A(7V+1) 

ae = e e a) = e ^ ^ ' 


a. 


Dirac’s operators 6^ are defined in terms of the aforesaid amplitude and angle variables: 

br = e--^®^y/Wr, [Qs, Nr] = ih5rs , 

whence one finds, using the representation 0^ = and consequently es®’’ = e 

decrease Nr by unity, whereas e will increase Nr by unity. 

a _ 

Thus, setting a = in (12.9p . we find: 

f{Nr) = f{Nr + or f{Nr) = f{Nr + l)es^’' . 


This brings us to the following set of useful formulae: 

br = y/Wr = x/Nr + Nr + le^ , 

bl = = e^^^yjNr + 1 = = e~^yjNr + 1, 

Nr = blhr . 


( 2 . 8 ) 

etc. 

(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 

d 

dNr will 


( 2 . 12 ) 


(2.13) 
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with Nr equal to zero or a positive integer for the number of particles in the state r. 

Let’s pick up the story at (12.5p and add to qn{t) the dimensionless conjugate momentum operator 


Pn{t) = 

V 2i 


(2.14) 


Then the Born-Jordan commutation relations in dimensionless form are given by 

[o-nj^rn] ^nm j ;Pm(^)] 

The transition from the quantum matter-field to the total energy operator of the string is given by 


H=ydx 

0 


a 


dip 

'm 


+ T 


dip 

dx 


^ ^ f -|- OnClji 


n=l 


oo . ^ 

^ ^ “t” “ 


71=1 


(2.16) 


With the following replacements of the operators by c-numbers 




V2?l/ 


1/2 


Cr)C 


we recover the classical expression (jl.l4l) for the total energy of the string. The real classical 
dynamical variables take the form 


h 


OJr. 


h 


UJr 


1/2 


quit) qn(t) = Cn COs(uJnt + 4>n) , 


1/2 


Pnit) Pn{t) = -UJnCn Sm{uJnt + 4>n) = (fnit) ■ 


:.cl { 


Now we are in a position to present a short introduction to the calculation of the mean-square 
energy fluctuation for the quantized string field. 

Let us first calculate the energy of a piece of length a of the string in the narrow frequency 
range a;, a; -|- Aw: 


m-\l 


dx 


dr\‘ (dr 


One obtains with H from (j2.16p : 

-^a,A;(^) — ^ ^ ^ \_QnQmI^nm PnP7n^nm\ 


n^m 


\—^ 1-/2 7 / \ \ 

+ 2^ fujJn [qn - Pn) sm(2A:na)— , Wn = ckn = y — jn, An = — . 


1 

n 

, _ sm{kn - km)a sm{kn + km)a 

“ kn- km kn + km 

sm{kn - km)a sin(/c„ -|- km)a 




kn kyyi 


kn H“ km 


(2.17) 


(2.18) 


11 
















Since we will be interested in the fluctuations in a narrow spectral range with Xn = ^ much smaller 
than both I and a, the final term in (I2.17D can be neglected due to fast oscillations, so the the 
deviation of the energy from the diagonal part, jH in (I2.17I] . is given by 




21 


^ ^ ^\/ \QnQm^nm PnPm^n 


(2.19) 


n^m 


From here on the reader can follow the explicit calculations outlined in chapter 5.2 of ref. 


• &• 


There 


is also Jordan’s original calculation in the Drei-Manner-Arbeit of Born, Heisenberg and Jordan 


B- 


Furthermore, it is rewarding to read up on the string field history in Heisenberg’s Chicago lecture 
of 1929 B Jordan’s conference contribution in Charkow in May of 1929 [s]- 
The result of Jordan’s work is 


AE2 = 




( 2 . 20 ) 


in complete analogy with Einstein’s mean square energy fluctuation formula for black body radia¬ 
tion: 






( 2 . 21 ) 


In either case, the origin of the particle-like fluctuation, second part in (I2.20p . is the presence of 
the zero-point energy term due to the non-commutativity of the quantized amplitudes qn{t) and 
Pm{i)- The appearance of this particle-like fluctuation is a kinematic effect which arises from the 
quantized nature of the string (photon) field. Jordan himself considered his fluctuation calculation 
as “almost the most important contribution I made to quantum mechanics.” In fact, it was the 
beginning of non-relativistic quantum field theory, if only for free fields. However, very shortly 
thereafter Dirac started his first attempt to formulate quantum electrodynamics. Nevertheless, 
Jordan introduced for the first time the quantum field operator (second quantization) as equivalent 
description to many particle wave functions. This is until the present day the accepted language 
for treating relativistic particle processes - with all its inherent problems that stem from products 
of relativistic local quantum field operators. 


III. DIRAC’S PAPER ON “THE QUANTUM THEORY OF THE EMISSION AND AB¬ 
SORPTION OF RADIATION” 


Jordan’s contribution in the final section of the Drei-Manner-Arbeit from November 15, 1925, 


and Dirac’s paper from February 2, 1927, [j] contain the foundation of quantum field theory, in 


12 
















particular of quantum electrodynamics. However, it is not correct to say that Dirac also invented 
the so-called second quantization procedure as we understand it today. Second quantization in its 
modern meaning as treating effectively many particle problems with the aid of quantum field opera¬ 
tors is Jordan’s contribution. Dirac baptized his procedure in this way because he introduced a pair 
of canonical variables, br and bl, which satisfy commutation rules and then act on a many-particle 
Schrodinger wave function, while Jordan elevated the one-particle Schrodinger wave function to 
charged matter-field operators of ordinary space-time. As “second quantized” field 

operators the ip{r,t) and satisfy certain equal-time commutation relations. We have al¬ 

ready seen in chapter 2 how these follow from the quantum action principle. It is only necessary to 
change the generator G(t) of the string theory to the appropriate form for the Schrodinger theory. 
Without further proof, the momentum conjugate to is given by n(r,t) = so that 

G{t) = j (f'r'ih'ip\r\t)5ilj{r't) , (3.1) 

which yields 

5'ijj{r,t) = J (3.2) 

From here we obtain 

= 5{r-r'). 

We also add the relations 

[V>(r,t),V’(r',f)] =0, 

= (3.3) 

Expanding and 'ijj^ in terms of a system of orthogonal c-number functions ifrif), 

= '^brit)ipr{r) , = '^bl{t)(p;{r) , j £’ripr{r)ipl{r) = 5r,s , (3.4) 

r u 

the commutation relations (|3.3p take the form 

bfb^g b\br — 5rs ) 

bj^bg bgbj' — 0 , 

blbl — blbl = 0 . (3.5) 

These equations are satisfied by (cf. also (|2.11l) and (|2.121) 1 

br = , bl = , (3.6) 
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where Nr and 0^ are Hermitian operators which satisfy the relation 

= f{Nr + l)e-H0’- . (3.7) 

The eigenvalue of the operator Nr are the integers = 0,1, 2,.... 

Nr = blbr is the number of particles in the r — th state. What is essentially new in the Dirac paper 
is that he considers the number Nr of systems in the r — th energy level as a canonical variable. In 
the space of such variables Dirac establishes the Schrodinger wave equation, i.e., the Schrodinger 
function becomes a function of the many Nr and time t. Let us see how it works. Consider an 
atom interacting with a radiation field. The Hamiltonian is given by 

H = Ho + V, (3.8) 


where H is the Hamiltonian of the unperturbed atomic system which is under the influence of the 
perturbation V. The wave function of the total system satisfies the time-dependent Schrodinger 
equation 


We then project onto the basis-states of the unperturbed system Hq, which are labeled by r with 
eigenvalues Wr 


• • , t\r)o = (•••; t\Ho\r)o = Wr{- ■ ■ ; f|r)o 

or ih-^ = iLoV'r = Writer ■ 

Ot 


(3.9) 


Expanding the Schrodinger wave function (• • • ;f|) in terms of (• • • ;t\r)o we have 


(••• ;i|) = ;ik)oo(?’l) 

r 

or = ''^^ar{t)ipr ■ (3.10) 

r 

ar{t) = (• • • ',t\r)o is the transformation function. \ar\‘^ is the probability of the system being in 
the state r at time t. The probabilities are normalized so that ^ = 1. To apply the theory 

r 

directly to an assembly of N similar independent systems, we must multiply each of these Ur hy^/N 
so as to make ^ |a,.p = A^. As a result we obtain in \ar\‘^ the probable number of systems in the 

r 

state r. Note that must be an integer. From (j3.9h and (I3.10h we obtain the rate of change of 
the expansion coefficients 

= '^Vrsds (3.11) 

S 
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and that for the complex conjugate a* 

= (3-12) 

s 

If we regard and iha* as canonical variables, the equations (13.111) and (I3.12|) take the Hamiltonian 
form with the Hamiltonian function 


Fr = '^<Vr, 


(3.13) 


namely, 


dur 1 dFi .f^da* dFi 


dt ih da* ’ 


dt dar 


(3.14) 


In order to express the complex canonical variables ar and a* in their real and imaginary parts, 
Dirac then introduced the real canonical variables Nr,4>r by the transformation 


= a* = _ 


In terms of these variables the Hamiltonian Fi becomes 


(3.15) 




(3.16) 


and the equations that determine the rate at which transitions occur have the canonical form 

dFi ■ dFi 


'd(l)r ’ 


dNr ■ 

Actually, Dirac found it more convenient to deal with the canonical variables 


(3.17) 


br = are-i^^\ K = 


(3.18) 


Wr being the energy of the stationary state r. We have equal to jarP, the probable number 
of systems in the state r. For br we find 


ihbr = Hrsbs = Wrbr + Vrsbs 


(3.19) 


where we put 


V -V 


— Wrdrs F Vrs • 


(3.20) 
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Taking the 6’s canonical variables instead of the a's means dealing with the Hamiltonian 


F = ^b;Hrsbs, 


(3.21) 


where we made a contact transformation 


b^ = ivysg-ie. 5* ^ ATi/yy- 


(3.22) 


to the new canonical variables Nr and 0^. The Hamiltonian (I3.21h will now become 


F = Y^ HrsN^^Ny^ 


e« 


(0r-0s) 


r,x 


The equations for the rate of change of Nr and 0^- will take the cononical form 


Nr = - 


dF 

dOr 


Qr = 


dF 

'm'r 


The Hamiltonian may be written 


(3.23) 


(3.24) 


F = Y, WrNr + Y, VrsN^N^ei 


{Br-Ba) 


(3.25) 


The first term, ^ WrNr, is the total proper energy of the assembly, and the second may be regarded 

r 

as the additional energy due to perturbation. 

Thus far Dirac had used only c-numbers. Now we assume the variables br,ihb* to be canonical 
g'-numbers satisfying the quantum condition 

brihbl — ihblbr = ih 


or 


brbl — b\.br = 1 


(3.26) 


and 


brbs - bgbr = 0 , blbl - = 0 , 

brbl — blbr = 0 (s y r). 


(3.27) 


The previous transformation equations must now be written in the form 


— -B 


br = e“S®’-iVy2 = {Nr + l)^/y 

bl = (ivyi/y:!®’- = e^®^{Nr + 1)^/^ 


(3.28) 

(3.29) 
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in order that Nr, &r be canonical variables satisfying the commutation rules 


[0r, Nr] = ih 


or, equivalently, 




-Pi 


(3.30) 


The equations (j3.26l 13.27113.28113.29P are precisely the same equations as the ones that are listed 
in (13.5113.6113.7p . The Hamiltonian becomes 


F = Y,blHrshs = Y,{Nrf/\i®^Hrs{Ns + 

r,s r,5 

= ^ ^ Vrs{Nrf^‘^{Ns + 1 - (3 31 ) 

r r,s 

and the Schrodinger equation can be written in terms of N^ by noting that and play 

the role of creation and annihilation operators. With the Schrodinger amplitude for observing 
particles at A^{, N 2 ,... N ^,... 


{N[,N', ... iv;,.. = HN[,N! 2 , ...N',...,t) 

we obtain 

,. ..N:.,...,t) = ^iN[,...,N' + l,...-t) 

,... iV;,..., t) = V;(iV(,..., fV; - 1,...; t). 


(3.32) 

(3.33) 


The evolution of ■ ■ ■ > ^) can thus be written as 

ih^^i;iN[,N',...;t) 

= FV>(iV(,iV',...;t) 

= ^W'N'i;iN[,N',...;t) 

+ ^u„(iV;)V2(iv' + 1 _ ... ,Ar; - 1,... ,iV' + 1,... ;t). (3.34) 

r,s 

Observe that so far we have a theory of free photons. This is because in the Lagrangian (13.211) we 
see that each creation operator bg is accompanied by an annihilation operator bl; hence there is a 
continuous transition of photons from one mode to another. Hence there is no hope of calculating 
the probabilities for emission or absorption of photons by an atomic electron. Dirac cures this 
situation by appending an extra interaction term to the original Hamiltonian. It takes the form 

Fin = ^ (^Vr4 + • 

r^O 


17 















With the use of such a linear interaction term, Dirac is able to compute Einstein’s famous A and 
B coefficients after quantizing the {b*,br) according to the {Nr,&r) variables: 

Hin = ^2 \/Nr + 1 + • 

r 

We will refrain from further elaboration of Dirac’s second quantization prescription simply be¬ 
cause his way of quantizing a field theory to become a theory of many particles which arise from 
excitations of underlying quantum fields is of mere historical interest now. The same fate was 
suffered by Dirac’s hole theory, which was replaced by a highly successful theory of electrons and 
positrons, quantum electrodynamics. It is here that the product of local second-quantized field 
operators, Jordan’s invention, proved to be the correct instrument to handle with great success all 
the interaction processes where electrons, positrons and photons are involved. 


IV. EQUIVALENCE OF JORDAN’S AND DIRAC’S SECOND QUANTIZATION PRO¬ 
CEDURE 

Let’s look at the relation between the usual one-particle Q.M. and the second quantization formal¬ 
ism. The connection is made at the level of the expectation value of a physical property /: 

b',b" 

(6"|) is the wave function of the state |) described with respect to the measurements h”. {\b') is the 
complex conjugate wave function described with respect to the measurements b'. From (14.ip we 
obtain 

= (4.2) 

b',b" 

This looks the same as the operator 

F=Y,a{h')HhV\b'')a{b"). (4.3) 

b',b" 

Here a{b')'^ and a{b") are operators. It is as though we had taken the numerical wave functions for 
a single system and extended them into operators, ending up with a description of a many-particle 
system. This is the origin of the name “second quantization”. To verify that our replacements 

{\b') ^ a{h')^ (4.4) 
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do not correspond to anything new but are merely a convenient description of n-particle systems 
(at the moment n = 1), let’s consider an arbitrary linear operator X acting on a one-particle state: 


X = \h'){h'\X\b"){h"\ ,i.e., X in5^ representation . 


(4.5) 


Then 


X\)oid = '^\b'){b'\X\b''){b''\)oid, (6"|)o«dwave function '(/’(&" 

b',b" 

The one-particle state can be described in either way: 




(4.6) 


or 


\)ne^ = Y.a{b"')^m{b"') , \h') = a(6')t|0). 


b'" 


Hence we obtain 


(4.7) 


X\)nen, = 16')(6'|X|6")(6"| a(6"')^|0)V’(6"') 

b'" 

= E |(>'>((>'|.>6|(>">(0|o((>")»('>"')'|0)?*((>"') 

b',b",b"' 

a{b”)a{b”'y = + a{b'")^a{b") 

a(6")|0) =0,(0|0) = 1 

= \h'){h'\x\b"mh") = Y, \h'){b'\x\h"){h"\)oid 

b',b'' b',b'' 

= X\)m. 

Now let the first particle have value bi, the secon 62, etc. and let every mode contain only one 
particle. This iV-particle basis state is then described by 

0(61)^0(62)^ ••• a(ftAr)'^IO). (4.8) 

Generalizing (14.71) to a general iV-particle state \(p') in which the particles have a wave function 
<y9(6i, 62,.. ., b]\f) we simply write a linear combination of states (14.8^ : 

\(j)') = Y ^{bi,b2,-■ ■ ,bN)a^{bi)a^b2).. .a^ibN)\0). (4.9) 
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If ip is properly symmetrized, then the amplitude for observing one particle in mode one in 
mode 62 ) etc. is given by 


{b[, b' 2 ,..., b'j^\(j)') not yet symmetrized. 

We want to check that the second quantized version of Jordan’s gives the same result as the one 
obtained by our old Schrodinger formulation. To show this, take another state 

\^') = 'fib[,b2,---,b'j^)a\b[)...a\b'N)\0) (4.10) 

and calculate the inner product 

i^'W) = X] <P*{b[,...,b'N)ip{bi,...,bN) 

6 ',... 6 '^ 

(0|o(6)v) • • • a{b[)a^bi)... ot(6Ar)|0). (4.11) 

Now we are going to apply the quantization relations (C.R.’s) [a{b'),a{b)^ = dyb to bring the a's 
to the right where they can act on | 0 ), a| 0 ) = 0 . 

To evaluate the vauum expectation value in (j4.11jl . let’s forget the problem for the time being and 
consider the function 

Y,a(b)a^h) 

e ,Oi{b) parametric function. (4.12) 


Now take this expression and differentiate: 


d d d 

da{bi) da{h 2 ) ''' da{bi<i) 


J2a(b)a'f(b) 

3 b 


a=0 


a'^(bi)a'^(b2}... a^bjy) ■ 


In the same way 


d 


d 


d 


da'{b\) da'{b'2) da'{b'j^) 


T.o^'ib')a{b') 

3 b' 


= a{b'^)a{b'2) ...a{b'j^) 


a '=0 


Then 


d 


d 


d 


d 


im- E i' (t'l, -X) *=('>!,..■ 


•( 0 |e 


6 i,... 6 jv 
6 ',... 6 '^ 

E «'(&')«(&') Y.<b)aHb) 


b' 


| 0 )q=q '=0 
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To compute the quantity (0| • • • |0) we write 


j:a'{b')a(b')j:a{b)at{b) -j:a'{b')a{b') j: a'{b')a{b') 

(Ole^-' e »' e>>' | 0 ) 

SF{a^)S-^ = F{Sa^S-^) 

just write in power series: S{a^)^S~^ = Sa^S~^S... Sa^S~ 

= (Sa^S-^)^ 

T,oi{b')a{b') -T,F{b')a(b') 

= (0|e*’' ^a(6)a'^(6)e >>' |0) . 

b 


To proceed, we have to know 


j:Fib')a{b') 

e*’' a^b)e>>' =F{a}. 


Taking the derivative of this expression, we first set up a differential equation: 


d ( Y.F{b')a{b') -j:F{b')a{b') 

‘ a\b)e 0' 


da'{b") 

Y.»'{V)a{V) , . . N -Y.F{b')a{b') 

= e*>' ^o(6'')o^(6) — a^(6)a(6")j e 


=<5: 


bb" 


= hb" = 


d 


da'{b"' 

From the last equation we obtain 


■(«'(&))• 


E«'(fe')<^(fe') , -E «'(&')«(&') 

F{a ) = e^' a*{b)e = a. {b) + const. 


and since F{a'{b')'\ — a'(b) is independent of a', we can determine its value at a' = 0 
with the constant of integration (6). From (I4.14p follows 

T.F{b')a(b') -j:F(b')a(b') 

e>>' a'{b)e *>' =a'(6) + o'(6) 


and from (14.131) : 


E«(fe)(<^t(&)+«'( 6 )), ^ E^W^hb) 

(0 |e'> |0 )=e‘- 


Thus, the scalar product looks like 


= X] <f>*{b'i,---,b'j^)ip{bi,...,bN) 


d 


d E«W«'(&) 


da'{b[) da{bN) 


b',...,b' 


a=a'=0 


(4.13) 


(4.14) 


to leave us 

(4.15) 


(4.16) 
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When the first N derivatives with respect to a{bi) are performed, we are left with 


(/9 symmetric function = A^! E ip*{bi,...,bN)^{bi,...,bN), 

bi,b2,---bM 

which, after redefinition of \(p): 

l</>) = V’{bi,---,bN)a\bi)...a\bN)\0) (4.17) 

yields just the old version of the scalar product. Now the amplitude for observing particles at 
b'i,b' 2 ,... ,b'j^ is given by 

(6i,...,6'jvl</>) = ;^^(0|a(6'i)...a(6'jv)l</>) = • • •, &jv) (4.18) 

m) = 1 (4.19) 

!</>)= X] \bu-■ ■ ,bN){bi,... ,bN\(l)) ■ (4.20) 

In other words, the operator 

1= ^ |6i,... ,67 v)(&i, ..., &Af| (4.21) 

fel,...,fejv 

is the unit operator when operating on properly symmetrized Wparticle states. 

Finally we come to time-dependent problems. Given H, then the equation of motion for the 
operator X{t) is given by {h = 1) 

^X{t) = -[X{t),H] 
at i 

with the solution 


X{t) = . 

f02 

If the single particle Hamiltonian is given by i7 = , then in a representation where p is diagonal, 

we obtain 


H = Y^a){p)^a{p) 

^a{p,t) = ^[a{p,t),H{t)] = -i^a{p,t) 


_ P^, 

a{p,t) = e a{p), 


a{p) = a{p, t = 0). 
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Likewise 


a){p,t) = e*^*a^(p). 

In general, for a 2nd quantized operator X, 

X = ^ a\p){p\X\p')a{p '), 

PV 

we have the time dependence 


which satisfies 


X(t) = ^ a\p,t){p\X\p')a{p',t) 
PP' 


-\x m 


For the operator H = ^ we have the equivalent expressions: 

X{t) = ^ a’^(p)e*^*(p|X|p')e“*^*a(p') 

p,p' 

= a^p){p\e^^^Xe~^^\p')a{p') 

PP' 

PP' 

Now we want to switch to the coordinate basis. Here we write instead of a{p) the operator 
a{x) = '0(®)- 


Since 


we obtain 


-iX p 


\x) = Y \p)(p\^) = Y 

p p V L 


and 


Therefore 


V^t(*) = ^at(p)^e-®'f 

= Y • 


{x),i)\y)\ = Y ^Hp),aHp')y^'Pe 


-iyp' 


PP' 

Y F —)• oo 


P 


d?p 

(27r)^ 


e^P<^-y) = 6\x - y) 
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Now the operator equation of motion is given by 


Substitute 


.d 




to obtain 




E 

P 

E 


\/F 


'^{x,t\p)a{p) 

P 




(Recall the Schrodinger equation for the amplitude {x,t\p)‘. 

d ^ 

= {xAH\p) = {x,t\^\p) 

X' 

But in the present chapter, denotes an operator. il>{x,t) destroys a particle at place x at time t. 


ACKNOWLEDGEMENT 

I enjoyed innumerable discussions with Nils Schopohl on the early days of quantum mechanics. 


[1] A. Einstein, Phys. Zeits. 10, 185 (1909) 

[2] M. Born, W. Heisenberg, and P. Jordan, Zeits. f. Phys. 35, 557 (1925) 

[3] S.S. Schweber, QED and the Men Who Made It, Princeton University Press 1994 

[4] P.A.M. Dirac, Proc. Roy. Soc. A114, 243 (1927) 

[5] P. Jordan and O. Klein, Zeits. f. Phys. 45, 751 (1927) 

[6] S. Varro, Fluctuation and Noise Letters Vol, 6, No. 3 (2006), World Scientific 

[7] W. Heisenberg, Physikalische Prinzipien der Quantentheorie, HI Taschenbcher/Mannheim (1958) 

[8] P. Jordan, Physik. Zeits. XXX, 700 (1929) 


24 



